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Abstract: In this paper optimal experimental designs for inverse quadratic regres- 
sion models are determined. We consider two different parameterizations of the 
model and investigate local optimal designs with respect to the c-, D- and E- 
criteria, which reflect various aspects of the precision of the maximum likelihood 
estimator for the parameters in inverse quadratic regression models. In particular 
it is demonstrated that for a sufficiently large design space geometric allocation 
rules are optimal with respect to many optimality criteria. Moreover, in numerous 
cases the designs with respect to the different criteria are supported at the same 
points. Finally, the efficiencies of different optimal designs with respect to various 
optimality criteria are studied, and the efficiency of some commorfly used designs 
are investigated. 
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systems, iJ— , c—, _D— optimality. 

1. Introduction 

Inverse polynomials define a flexible family of nonlinear regression models 
which are used to describe the relationship between a response, say Y , and a 
univariate predictor, say u [see eg. Nelder (1966)]. The model is defined by the 
expected response 

It 

where Pniu, 0) is a polynomial of degree n with coefficients 9q, . . . ,6n defining the 
shape of the curve. Nelder (1966) compared the properties of inverse and ordi- 
nary polynomial models for analyzing data. In contrast to ordinary polynomials 
inverse polynomial regression models are bounded and can be used to describe 
a saturation effect, in which case the response does not exceed a finite amount. 
Similarly, a toxic effect can be produced, in which case the response eventually 
falls to zero. 
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An important class of inverse polynomial models are defined by inverse 
quadratic regression models, which correspond to the case n = 2 in (II. ip . These 
models have numerous applications, in particular in chemistry and agriculture 
[see Ratkowski (1990), Sparrow (1979a, 1979b), Nelder (1960), Serchand, Mc- 
New, Kellogg and Johnson (1995) and Landete-Castillejos and Gallego (2000) 
among others]. For example. Sparrow (1979a, 1979b) analyzed data from several 
series of experiments designed to study the relationship between crop yield and 
fertilizer input. He concluded that among several competing models the inverse 
quadratic model produced the best fit to data obtained from yields of barley and 
grass crops. Similarly, Serchand et al. (1995) argued that inverse polynomials 
can produce a dramatically steep rise and might realistically describe lactation 
curves. 

While much attention has been paid to the construction of various optimal 
designs for the inverse linear or Michaelis Menten model [see Song and Wong 
(1998), Lopez-Fidalgo and Wong (2002), Dette, Melas and Pepelyshev (2003), 
Dette and Biedermann (2003), among many others], optimal designs for the in- 
verse quadratic regression model have not been studied in so much detail. Cobby, 
Chapman and Pike (1986) determined local D-optimal designs numerically and 
Haines (1992) provided some analytical results for D-optimal designs in the in- 
verse quadratic regression model. In particular, in these references it is demon- 
strated that geometric allocation rules are D-optimal. The present paper is 
devoted to a more systematic study of local optimal designs for inverse quadratic 
models. We consider the c-, -D-, Di- and £^-optimality criterion and determine 
local optimal designs for two different parameterizations of the inverse quadratic 
regression model. In Section 2 we introduce two parameterizations of the in- 
verse quadratic regression model and describe some basic facts of approximate 
design theory. In Section 3 we discuss several c-optimal designs. In particular 
Di-optimal designs are determined, which are of particular importance if dis- 
crimination between an inverse linear and inverse quadratic model is one of the 
interests of the experiment. As a further special case of the c-optimality crite- 
rion we determine optimal extrapolation designs. Section 4 deals with the local 
-D-optimality and £'-optimality criterion. It is shown that for all criteria under 
consideration geometric designs are local optimal, whenever the design space is 
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sufficiently large. We also determine the structure of the local optimal designs in 
the case of a bounded design space. These findings extend the observations made 
by Cobby, Chapman and Pike (1986) and Haines (1992) for the D-optimality cri- 
terion to other optimality criteria, different design spaces and a slightly different 
inverse quadratic regression model. 

2. Preliminaries 

We consider two parameterizations of the inverse quadratic regression model 

E{Y\u)=r,{u,9) , (2.1) 

where 6 = {Oq, 61,02)'^ denotes the vector of unknown parameters and the ex- 
pected response is given by 

u 

6*0 + 6li"u + (92^2 

or 



V2{n,e)= J'l ^ (2.2b) 

The explanatory variable varies in the interval lA = [s,t], where s > and 
< s < t < 00, or in the unbounded set U = [s, 00) with s > 0. The assumptions 
regarding the parameters vary with the different parameterizations and should 
assure, that the numerator in (j2.2ap and ()2.2bp is positive on lA. Under such 
assumptions the regression functions have no points of discontinuity. Moreover, 
both functions are strictly increasing to a maximum of size {9i + 2\/ 0q62)~^ at 
the point Umaxi = a/ ^0/^2 for parameterization ()2.2ap and to a maximum of 
size ^0(1 + V^i^)^^ at the point ?Xmax2 = \f^\[^ for parameterization (|2.2bp 
and then the functions are strictly decreasing to a zero asymptote. A sufficient 
condition for the positivity of the numerator is 6*0,^2 > 0, | | < 2\/^o^2 for 
model ¥12^ and 60,61,62 > 0, 2^/6^ > 1 for model ([23b]) . respectively. We 
assume that at each u £ U a normally distributed observation is available with 
mean rj{u,6) and variance a"^ > 0, where the function rj is either rji or 7J2, and 
different observations are assumed to be independent. An experimental design 
^ is a probability measure with finite support defined on the set U [see Kiefer 
(1974)]. The information matrix of an experimental design is defined by 



M{^,6)= I f{u,6)f{u,6)di{u), (2.3) 
u 
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where 

f{u,e) = ^r^{u,9) (2.4) 

denotes the gradient of the expected response with respect to the parameter 
9. For the two parameterizations (|2.2ap and (|2.2b|l the vectors of the partial 
derivatives are given by 

and 

f2{u,9) = ^^^^^^^^2 {^^-0^ + u + 92U^'~9i + u + 92U^) ' ^'^'^^ 
respectively. 

If N observations can be made and the design ^ concentrates mass Wi at the 
points Ui, i = l,...,r, the quantities WiN are rounded to integers such that 
X]^=i ^« — ^ [s^^ Pukelsheim and Rieder (1992)], and the experimenter takes 
rii observations at each point Ui, i = 1, . . . ,r. If the sample size N converges 
to infinity, then (under appropriate assumptions of regularity) the covariance 
matrix of the maximum likelihood estimator for the parameter 9 is approximately 
proportional to the matrix (a'^ /N)M~^{(^,9), provided that the inverse of the 
information matrix exists [see Jennrich (1969)]. An optimal experimental design 
maximizes or minimizes an appropriate functional of the information matrix or 
its inverse, and there are numerous optimality criteria which can be used to 
discriminate between competing designs [see Silvey (1980) or Pukelsheim (1993)]. 
In this paper we will investigate the -D-optimality criterion, which maximizes 
the determinant of the inverse of the information matrix with respect to the 
design ^, the c-optimality criterion, which minimizes the variance of the maximum 
likelihood estimate for the linear combination c^9 and the i?-optimality criterion, 
which maximizes the minimum eigenvalue of the information matrix M{^,9). 

3. Local c-optimal designs 

Recall that for a given vector c G M'^ a design is called c-optimal if 

the linear combination (F9 is estimable by the design that is Range(c) C 
Range(Af(^c5 ^))5 and the design minimizes 

c^M-{i,9)c (3.1) 
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among all designs for which c^9 is estimable, where M~{^,6) denotes a gener- 
alized inverse of the matrix M{^,6). It is shown in Pukelsheim (1993) that the 
expression (13. ip does not depend on the specific choice of the generalized inverse. 
Moreover, a design is c-optimal if and only if there exists a generalized inverse 
G of M(^c, &) such that the inequality 

{f'{u,e)Gcf <(^M-{i,,e)c (3.2) 

holds for all n G W [see Pukelsheim (1993)]. A further important tool to deter- 
mine c-optimal designs is the theory of Chebyshev systems, which will be briefly 
described here for the sake of completeness. 

Following Karlin and Studden (1966) a set of functions {go, . . . defined 
on the set U is called Chebychev-system, if every linear combination ^11=0 ^^iQii^) 
with X^ILo ^1 ^ *-* most n distinct roots on lA. This property is equivalent 

to the fact that 

det(5f(no), . . . , g{un)) / (3.3) 

holds for all uq, . . . ,Un ^ U with Ui 7^ uj {i 7^ j), where g{u) = (goiu), ■ ■ ■ , Qni'u))'^ 
denotes the vector of all functions [see Karlin and Studden (1966)]. If the func- 
tions go, . . . ,gn constitute a Chebyshev-system on the set U, then there exists a 
unique "polynomial" 

n 

(t>iu) ■.= ^a*gi{u) (aS,...,a: gM) (3.4) 

with the following properties 

(i) I (piu) I < 1 yueu 

(a) There exist n + 1 points sq < ■ • • < Sn such that (j){si) = (— 1)"""* 
for i = 0, . . . , n. 

The function (j){u) is called the Chebychev-polynomial, and the points sq, . . . ,Sn 
are called Chebychev-points, which are not necessarily unique. Kiefer and Wol- 
fowitz (1965) defined the set A* C as the set of all vectors c E 
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satisfying 



goiXn) Co 
gi{Xn) Cl 



(3.5) 



gn{xi) ... 9n{Xn) Cn 

whenever the points xi,...,Xn G U are distinct. They showed that for each 
c € A* the c-optimal design, which minimizes 

\ -1 

g{u)g'^{u)d^{u) 



U 

among all designs on lA, is supported by the entire set of the Chebychev-points 
sO) • • • , The corresponding optimal weights i^q, . . . , u)* can then easily be 
found using Lagrange multipliers and are given by 



En I I 



0, . . . ,n, 



(3.6) 



where the vector v is defined by 

V = {xx^y^xc, 

and the (n + 1) x (n + l)-matrix X is given by X = {gjisi))^j^Q [see also 
Pukelsheim and Torsney (1991)]. 

In the following discussion we will use these results to determine local op- 
timal design for two specific goals in the data analysis with inverse quadratic 
regression models: discrimination between inverse linear and quadratic models 
and extrapolation or prediction at a specific point Xe- We will begin with the 
discrimination problem, which has been extensively studied for ordinary poly- 
nomial regression models [see Stigler (1971), Studden (1982) or Dette (1995), 
among many others]. To our knowledge the problem of constructing designs for 
the discrimination between inverse rational models has not been studied in the 
literature. We consider the inverse quadratic regression model ()2.2al) and are 
interested in determining a design, which can be used to discriminate between 
this and the inverse linear regression model 



r]{u,e) 



9o + Oiu' 
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The decision, which model should be used could be based on the likelihood ratio 
test for the hypothesis Hq : ^2 = in the model (I2.2ap . and a standard calculation 
shows that the (asymptotic) power of this test is a decreasing function of the 
quantity (j3.ip . where the vector c is given by c = (0,0,1)-^. Thus a design 
maximizing the power of the likelihood ratio test for discriminating between the 
inverse linear and quadratic model is a local c-optimal design for the vector 
c = (0,0, 1)"^. Following Stigler (1971) we call this design local Di-optimal. Our 
first results determine the local Di-optimal design for the two parameterizations 
of the inverse quadratic regression model explicitly. 



Theorem 3.1 The local Di-optimal design for the inverse quadratic regres- 
sion model h2. 2a\) on the design space lA = [0, 00) is given by 



i §0 

Wq 



Wi 



1 — Wq — Wi 



(3.7) 



with weights 



Wo 
Wl 



92^0 + ^1' 



9op') /[(l + p)A] 



and A = 00(^0^2(1 + 6/>2 + p^) + 29ip{eip + V9^{l + pf)) ■ 
The geometric scaling factor p is defined by 

p = />(7) = 1 + ^ + ^2(1 + ^2) + (2 + ^2)7 + ^ (3.8) 

with 7 = 6i/\/0q62- This design is also local Di-optimal on the design space 
U = [s,t] (0 < s < t), if the inequalities < s < p^^ \/0q/92 and t > p^J 60/62 
are satisfied. 

The local Di-optimal design on the design space U = [s,t] for model 12. 2a\) is of 
the form 

t* _ ( ^ u'l ^2 
^■Di - , , . , , 

\Wq W\ 1 — Wq — W^ 
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if the inequalities s > p ^ cLnd t > p\J ^0/^2 hold, of the form 



Uq u'l t 



^k = [ : > , J' (3-10) 



if the inequalities s < p ^ y/Oo/62 and t < p^/ 60/62 are satisfied, and is of the 
form 

( s u'l' t \ , ^ 

Cd. = \ h 3.11 

' \w'^' w'i' 1 - < - w'^l 

if the inequalities s > p^^ \/6o/62 and t < p^ 60/62 hold. 
Proof: The proof is performed in three steps: 

(A) At first we identify a candidate for the local Di-optimal design on the 
interval [0, 00) using the theory of Chebyshev polynomials. 

(B) We use the properties of the Chebyshev polynomial ()3.4p to prove the local 
Di-optimality of this candidate. 

(C) We consider the case of a bounded design space and determine how the 
constraints interfere with the support points of the local optimal design on 
the unbounded design space. 

(A): Let /(u, 6) be the vector of the partial derivatives in parameterization ()2.2ap 
defined in (12. Sh . It is easy to see, that the components of the vector fi{u, 6), say 
{ fio{u,6), fii{u,6), fi2{u,6) }, constitute a Chebyshev-system on any bounded 
interval [s, t] C (0, 00). Furthermore for yo, yi > with yo / yi we get 



fioiyo,0) fio{yi,6) 
/ii(yo,^) fii{yi,o) 
/i2(yo,^) /i2(yi,^) 1 



and it follows that the vector (0,0,1)"^ is an element of the set A* defined in 
(13. 5p . Therefore we obtain from the results of Kiefer and Wolfowitz (1965), that 
the local Di-optimal design is supported on the entire set of Chebyshev-points 
{uq,u\.,U2} of the Chebyshev-system {fio{u,6), fii{u,6), fi2{u,6)}. If the sup- 
port points are given, say uq, ui,U2 the corresponding weights can be determined 
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by (j3.6p such that the function defined in (j3.1|) is maximal. 
Consequently the Di-optimality criterion can be expressed as a function of the 
points UQ,ui,U2, which will now be optimized analytically. For this purpose we 
obtain by a tedious computation 

^) '■= StT^ =^o(^o - - ^2M/N (3.12) 

where M(^, 9) denotes the matrix obtained from M(^, 9) by deleting the last row 
and column, u = {uo,ui,U2),9 = {9q,9i, ^2) and 

= (400^*0^^1 (^1 + 02Ui)u2 + 9q{ui{u2 - ui) + uo{ui + U2)) 

+ uouiU2{29lui + 29i92{uo{ui - U2) + ui{ui + U2)) 

+ ^li'^oiui - U2) + Uo{ui - U2)U2 + Ui{u\ + ul)))f. 

The support points of the Di-optimal design are obtained by maximizing the 
function T{u,9) with respect to uo,ui,U2. The necessary conditions for a maxi- 
mum yield the following system of nonlinear equations 
dT 

- — {u,9) =49oului{9i + 92Ui)u2 + 9'^{-2uoUi{ui - U2) 

OUq 

+ ul{ui - U2) + Uq{ui + U2)) + uluiU2{29jui + A9i92u\ 

+ 92{2uoUi{ui - U2) + Uq{u2 - ui) + u\{ui + U2))) • -Ri = 0, 

(3.13) 

dT 

- — {u, 9) = A9qUqu\u2[-{92u\) + 92U0U2 + 01 (no - 2ui + U2)) 

OUi 

+ 9q{uI{ui - U2)^ - 2uqUi{u\ + U1U2 - u\) + Uq{u\ + 2uiU2 - U2)) 
- Uqu\u2{29'1{u\ - UQU2) + A9i92Ui{uq{ui - 2U2) + U1U2) 

+ (^'ii'^oiui — U2)'^ + 2uoUi{ul — U1U2 — u\) 

+ ul{-ul + 2uiU2 + ul))) ■ R2 = Q, (3.14) 

dT 

- — {u, 9) = A9ququi{9i + 02^*1)^^2 + 9l{ui{ui - ^2)^ + uq{-u\ + 2uiU2 + ^2)) 

UU2 

+ UQUiul{29lui + 49i92uj + 6i|(uo(ui - U2f 

+ ui{ul + 2uiU2-ul)))- R3 = 0, (3.15) 

where Ri, R2 and R3 are rational functions, which do not vanish for all uq, ui, U2 
with < no < ui < U2. In order to solve this system of equations, we assume 
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that 

uo = — , U2 = r ■ ui (3.16) 
r 

holds for some factor r > 1, which wih be specified later. Inserting this expression 
in (13.14P provides as the only positive solution ul = ■s/QofO^. Substituting this 
term into ()3.13p or (j3.15p yields the following equation for the factor r 

201 (01 + 4y0^)r2 - 0002(1 - 4r - 2r'^ - ^r^ + r^) = 

with four roots given by 

ri/2 = 1 ± ± ^2(1 + ^/2) + (2 + 72)7 + ^, (3.17) 

r3/4 = 1 ± T y2(l + ^/2) + (2 + y2)7 + ^, 

where 7 = ^ ■ The factor r has to be strict greater 1 according to our 

assumption on the relation between no,ui and U2- This provides only the first 
solution in (j3.17p and the geometric scaling factor is given by (jS.Sh . Therefore 
it remains to justify assumption (j3.16p . which will be done in the second part of 
the proof. 

(B) Because the calculation of the support points p^^y^o/^i \/^o/^2 and 
P\l ^0/^2 in step (A) is based on assumption (j3.16p . we still have prove that these 
points are the support points of the local Di-optimal design. For this purpose 
we show that the unique oscillating polynomial defined by ()3.4p attends minima 
and maxima exactly in these support points. Recall that the vector of the partial 
derivatives of the regression function /i(n, 0) = {fiQ{u,6), fu{u,6), fi2{u,6)) is 
given by ()2.5p . We now define a polynomial t{u) by 

tiu) = /io(n, 0) + ai/ii(n, 0) + a2/i2(^^, 0) (3.18) 

and determine the factors ai and 02 such that it is equioscillating, i.e. 

t'{u*) = i = 0,1,2 (3.19a) 

t{u*) = c{-iy-'^ i = 0,1,2 (3.19b) 

for some constant c G M. By this choice the polynomial t{u) must be proportional 
to the polynomial 0(n) defined in (j3.4p . For the determination of the coefficients 
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we differentiate the polynomial t{u) and get 

-(6*0(1 + 2uai + 3ti^a2)) + ti(6'i(l - u'^a2) + 92u{3 + 2uai + u^a2)) 



t'{u) 



+ u{ei + e2u)f 



Substituting the support points u\ = -\/^o and 
obtain from (|3.19a|) two equations 



(3.20) 

1*2 = in ([3:20]) we 








72(72 



+ 



2{0iP + VW(3p2 _ 1) + 20op(p' - l)ai) 



V C'2C'2 

Wop'i-ieip) + VWiip" - 3))a2 



72^2 



The solution with respect to ai and 02 is given by 



29op 



which yields for the polynomial t{u) and its derivate 

n(-20op + + /0^)^i - pu{Oi + 2d2u)) 



tiu) 
t'{u) 



26op{eo + u{9i + 92u)y 



2PU) 



(3.21) 



0op{eo + u{ei + e2u)f 

respectively. A straightforward calculation shows that the third support point 



Uq = p~^\/0q/62 satisfies t'(Mo) = and that the three equations in (|3.19bp 
are satisfied. Therefore it only remains to prove, that the inequality | t{u) \ < c 
holds on the interval [0,oo). In this case the polynomial t{u) must be propor- 
tional to the equioscillating polynomial (j){u) and the design with support points 
p^^ \/0q/62, \/Oq/02 and py/Oo/92 and optimal weights is local Di-optimal. 
Observing the representation p.2ip shows that the equation t'{u) = is equiva- 
lent to 



with roots 



02u)WOoP - V92u){Vdo + V92u){Vdo - V92Pu) = (3.22) 



no 



1/^0 / ^0 J 9o 
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Therefore the function t{u) has exactly three extrema on R+. Furthermore if 
ti ^ oo, we have t{u) — > and it fohows that | t{u) \ < c holds for all n > 0. 
Consequently, the functions t{u) and (j){u) are proportional and the points Uq = 
p^^ y/Oo/62, u* = y/Oo/62, U2 = P\/ 00 /O2 are the support points of the local Di- 
optimal design. The explicit construction of the weights wq and wi is obtained 
by substituting the support points Uq, ul and U2 into ()3.6p . 

(C) We finally consider the cases (13.9p . (|3.10p and (|3.1ip in the second part 
of Theorem 13. 1^ which correspond to a bounded design space. For the sake of 
brevity we restrict ourselves to the case (|3.9p , all other cases are treated similarly. 
Obviously the assertion follows from the existence of a point Uq > 0, such that the 
function T{uo,ul,U2,0) is increasing in uq on the interval (0,Uq) and decreasing 
on {u*,ul). 

For a proof of this property we fix ui,U2, and note that the function T{uq) := 
T{uq, ui,U2, 0) has minima in uq = and uq = ui, since the inequality T{uq) > 
holds for all uq G [0, ui] and T(0) = T[ui) = 0. Because T[uq) is not constant, 
there is at least one maximum in the interval (0, ui). In order to prove that there 
is exactly one maximum, we calculate 

T'{uq) = ^{uq, ui,U2, 6) = 2uo{uo - ni)uf (m - U2ful ^^'^'^\ , (3.23) 
duo Pdiuo) 

where Pg is a polynomial of degree 9 (which is in the following discussion without 

interest) and the polynomial P4 in the numerator is given by 

P4(^io) = 46*0^0^1 (6*1 + 62Ui)u2 + 6l{-2uoUi{ui - U2) + uf{ui - U2) 

+ u^{ui + U2)) + UqUiU2{29'Iui + A6i92u{ + 92{2uqui{ui - U2) 

+ UI{U2 - Ul) + u\{ui + U2))). 

The roots of the function T" are given by the roots of the polynomial P4. Differ- 
entiating this polynomial yields the function 

(Uo) = WQUQUi{di + 92Ui)U2 + 29l{ui{u2 - Ul) + Uq{ui + U2)) 

OUq 

+ 2uoUiU2{29lui + A9i92ul + 9l{-2ul{ui - U2) 

+ 2>uoui{ui - U2) + u\{ui + U2))), 
which has only one real root. Consequently ^4(^0) has just one extremum and 
therefore at most two roots. The case of no roots has been excluded above. If 
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P4,{uo) would have two roots, then the function T{uo) has at most two extrema 
in the interval {0,ui). However, the function T{uq) is equal to zero in the two 
points and ui and in the interval (0,-ui) strictly positive. Therefore the num- 
ber of its extrema has to be odd and T(uq) has exactly one maximum on (0, ui), 
which is attained for given (ni,M2) = iul,U2) at a point Uq G (0, mJ). 
Assume that the design space is of the form U = [s,t]. If the inequality s < Uq 
holds, ()3.7p remains the local Di-optimal design. However if the inequality s > Uq 
holds, the function T[uq) is maximal in s, and it follows that (j3.9p is the local 
Di-optimal design. 

Remark 3.1 Note that part (A) of the proof essentially follows the arguments 
presented in Haines (1992) for the D-optimality criterion, who considered the 
model 

Po + /?!(« + «)+ f32{u + ay 
However, the proof presented by Haines (1992) is not complete, because she did 
neither justify the use of the geometric design, nor proves that the system of 
necessary conditions has only one solution. In this paper we present a tool for 
closing this gap, as demonstrated in part (B) of the preceding proof. 
It is also worthwhile to mention that an analogue of Theorem 13. II does not hold in 
the four-parameter model discussed in Haines (1992). For example if /3o = /32 = 
/?! = —1.8 and a = 0.1 we obtain by numerical computation that the local Di- 
optimal design is supported at the Chebyshev-points { 0,0.6272,0.9861, 1.8714 } 
and there does not exist a similar geometric spacing behaviour as in the models 
considered in this paper. □ 

The following theorem states the corresponding results for the inverse quadratic 
regression model with parameterization (|2.2bj) . The proof is similar to the proof 
of the previous theorem and therefore omitted. 

Theorem 3.2 The local Di-optimal design for the inverse quadratic regres- 
sion model ^2. 2h\) on the design space U = [0, oo) is given by 
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with 

WO =(7^01 + y^ip + y%eip^f{l + + P))/[{1 + P)A] 

w;i =(201 + + V^(l + P'))/A 

A = 0i(p(2p + 3v/^(l + p)^) + 01^2(1 + 2v^(l + p)\l + P^) 
+ p(8 + p(6 + /9(8 + p))))). 

r/ie geometric scaling factor p is given by 113. 8\) with 7 = {V0i92)^^ ■ This design 
is also local Di-optimal on the design space U = [s,t] (0 < s < t), if the inequal- 
ities < s < p^^ \fQ\fO2 O'lT'd t > p\J 61/62 are satisfied. 

The local Di-optimal design on the design space lA = [s, t] for the inverse quadratic 
regression model 112. 2b\) is of the form i3. 9\) if the inequalities s > p~^ \f6\l62 
and t > p-y/ 61/62 hold, of the form ^3.10\} if the inequalities s < p^^ \J 61/62 
and t < p\j6\/62 are satisfied and is of the form i3.11]) if the inequalities 
s > p~^^J 61/62 and t < p^ 6^/62 hold. 

In the following discussion we concentrate on the problem of extrapolation 
in the inverse quadratic regression model. An optimal design for this purpose 
minimizes the variance of the estimate of the expected response at a point and 
is therefore c-optimal for the vector Cg = fi{xe, 6) in the case of parameterization 
()2.2ap . and for the vector Cg = f2{xe,6) in the case of parameterization (I2.2b|) . 
respectively. If Xe is an element of the design space U it is obviously optimal to 
take all observations at the point Xe, and therefore we assume for the remaining 
part of this section that U = [s,t], where < s < t and < x^ < s or Xe > t. 
The following result specifies local optimal extrapolation designs for the inverse 
quadratic regression model which are called local Cg-designs in the following dis- 
cussion. The proofs are similar to the proofs for Di-optimality and therefore 
omitted. 

Theorem 3.3 Assume that U = [s,t], where < s < t and < x^ < s or 
Xe > t, and let p denote the geometric scaling factor defined in i3. 8]) with 7 = 
6i{\/6q62)~^ . If < s < p~^\/6o/62 and t > py^ 60/62, then the local Ce-optimal 



OPTIMAL DESIGNS FOR INVERSE REGRESSION MODELS 

design ^* for the inverse quadratic regression model 112. 2a\) is given by 
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1 [K IK ^ 

Wo Wl 1 — Wq — Wl 



(3.25) 



where 



Wo =(7^-Xex/^)(-Xex/^+ y^p) (Ooy^ + 6iy%p + 9oy^p' 

x((l + p)A)-i 

Wl =(29oyT2 + ^1 V^)'p( - XeV^ + y%p) - XeV0'2p)/X 



with 

X = 00 (0^^2(1 + + /) + 60(26^ + 2eip(v^(l + p)' - 4Xe02(l + p')) 

+ 02Xe( - 2v^(l + p)'(l + p') + Xe02(l + Op^ + p^))) 

+ 0iXep(2y^e2Xe(l + p)' " + p( " 2x^02 + \/^(2 + p))))) 

T/ie local c^-optimal design for the inverse quadratic model 12. 2a\) is of the form 
113. 9\) if the inequalities s > \J Q0IO2 o-nd t > Pa/^o/^2 hold, of the form h3. 1 0\) 
if the inequalities s < p^^ \fQ0IQ2 o-nd t < py^o/^ are satisfied and of the form 
113. if the inequalities s > p^^\/Oo/92 and t < p^ O0/O2 hold. 



Theorem 3.4 Assume that lA = [s,t], where < s < t and < Xg < s or 

Xg > t, and let p denote the geometric scaling factor p defined in /i3. 8\) with 
7 = {^/6^)-\ IfO<s< p"V^i/^2 and t > py^6i/92, then the local Cg- 
optimal design for the inverse quadratic regression model i2. 2b]) on the design 
space U = [0, 00) is given by 




(3.26) 
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el[el92{l + 6p^ + p^)+Xep(^-^/0^2 + 2^/0^2O2Xe{l + p)'' 

-p{- 2Xe7^+ 7^(2 + />))) + 0i(2p2 + 2p{y^{l + p)^ 
- 4XeV^(l + + Xe{ - 27^^2(1 + + P^) 

+ xA'(l + 6/,2^p^)))) 

// t/ie design space is given by a finite interval [s,t], < s < t, then the local 
Cg-optimal design for model h2. 2a\) is of the form h3.9^) . if the inequalities s > 
p^^ \/0i/62 and t > p-\fW\jQ2 hold, of the form Ii3.10\) . if the inequalities s < 
p^^\f6\Q2 andt < p^J 61/62 are satisfied, and of the form Ii3.11\) if the inequalities 
s > p"V^i/^2 and t < p^ 6x1 62 hold. 

Note that for a sufficiently large design interval all designs presented in this 
section are supported at the same points, the Chebyshev points corresponding 
to the Chebyshev system of the components of the gradient of the regression 
function. In the next section we will demonstrate that these points are also the 
support points of the local E'-optimal design for the inverse quadratic regression 
model. 

4. Local D- and i?-optimal designs We begin stating the corresponding 
result for the D-optimality criterion. The proof is omitted because it requires 
arguments which are similar as those presented in Haines (1992) and in the proof 
of Theorem 13.11 



with 

X = 



Theorem 4.1 The local D- optimal design for the inverse quadratic regression 
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model i2. 2a\) on the design space lA = [0, cxd) is given by 



eD= I I ' ^^-^^ 

3 3 3 



with the geometric scaling factor 

P= 2 ' 



(4.2) 



where the constants 6 and 7 are defined by 5 = (1/2) (7 + 1 + y^-f^ + 67 + 33) 
and 7 = Ol{^/9o92)^^ , respectively. This design is also local D-optimal on the 
design space U = [s,t] (0 < s < t), if the inequalities < s < p^^ ^/Oq/Oq, and 
t > p\J 6^/62 are satisfied. 

The local D-optimal design on the design space lA = [s, t] for the inverse quadratic 
regression model /(2. 2b\) is of the form i3. 9\) . if the inequalities s > p^^ \/0i/62 
and t > p\fW\J02 hold, of the form /i3.10\} . if the inequalities s < p^^ \/0i/92 
and t < p\/9i/92 are satisfied, and is of the form liS.ll]) . if the inequalities 
s > p~V^i/^2 and t < p^jQxlQ2 hold. 



Theorem 4.2 The local D-optimal design for the inverse quadratic regression 
model 112. 2b]) on the design space U = [0, 00) is given by 




with the geometric scaling factor p is given by (l§ with 7 = {V0^y\ This 
design is also D-optimal on the design space U = [s,t] (0 < s < t), if the 
inequalities < s < p~^ y/Oi/62 and t > py^ 61/62 are satisfied. 
The local D-optimal design on the design space U = [s, t] for the inverse quadratic 
regression model ^2. 2b\) is of the form liS. 9\) . if the inequalities s > p^^ \f6\l62 
and t > p-sj 61/62 hold, of the form \3.10\} . if the inequalities s < p^^ \J 61/62 
and t < p\/6i/62 are satisfied, and is of the form 113. if the inequalities 
s > p"V^i/^2 and t < p^j6x/62 hold. 



We will conclude this section with the discussion of the i?-optimaIity crite- 
rion. For this purpose recall that a design is local E'-optimal if and only if 
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there exists a matrix E G conv(S') such that the inequahty 

f'{u,e)Ef{u,9)<X^i^ (4.4) 

holds for sdl u G U, where Amin denotes the minimum eigenvalue of the matrix 
M{Ce,0) and 

S" = { zz' I ||2:||2 = 1, z is an eigenvector of M{^e,0) corresponding to Amin } • 

(4.5) 

The following two results specify the local £J-optimal designs for the inverse 
quadratic regression models with parameterization (j2.2ap and ()2.2bp . Because 
both statements are proved similarly, we restrict ourselves to a proof of the first 
theorem. 

Theorem 4.3 The local E-optimal design for the inverse quadratic regression 
model 112. 2a\) on the design space U = [0, co) is given by 



- Wo — Wi 




where the weights vuq, wi are given by i3. 6\) and c is the vector with components 
given by the coefficients of the Chebyshev polynomial, that is 

VO~o{291p^ + 2V^0i V^p(l + p)2 + 00^2(1 + 6p2 + 



elpji + p)^ + 8V9^eiV9^p{i + p2) + 2^0^2(1 + p)'(i + p^) 

(-l + p)2p 

V0^i20lp^ + 2^eiV0^p{l + p)2 + ^0^2(1 + + p^))^T 



V^(-l + p)V 

The geometric scaling factor is given by 113. 8\) with 7 = ^i(\/^o^)^^. This design 
is also local E-optimal on the design spaceU = [s,t] (0 < s < t), if the inequalities 
< s < p^^ \/0q/62 and t > py^o/^ are satisfied. 

The local E-optimal design on the design space lA = [s,t] for model 1^2. 2a\) is of 
the form ^3.9\) . if the inequalities s > p^^ \/0q/62 and t > py^ 60/62 hold, of the 
form i3.10\) . if the inequalities s < p~^ \/6o/62 o-nd t < p^ 60/62 are satisfied, 
and of the form \3.11\) . if the inequalities s > p~^ \/6o/62 and t < py^ 60/62 hold. 
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Proof: It is straightforward to show that every subset of {fioiu, 9), fu{u, 6), 
fi2{u,6)}, the components of the vector fi{u,9), which consists of 2 elements, 
is a (weak) Chebychev-system. Therefore it fohows from Theorem 2.1 in Imhof, 
Studden (2001) that the local i?-optimal is supported at the Chebyshev points. 
The assertion regarding the weights finally follows from (j3.6p observing that the 
results of Imhof and Studden (2001) imply that the local E'-optimal design is 
also c-optimal for the vector c with components given by the coefficients of the 
Chebyshev polynomial. □ 

Theorem 4.4 The local E -optimal design for the inverse quadratic regression 



model \2. 2h\) on the design space lA = [0, oo) is given by 




where the weights wq, wi are given by i3. 6\) and c is the vector with components 
given by the coefficients of the Chebyshev polynomial, that is 

. 2(2p + + pfKp + vma + p^)) 

(-1 +/9)V 

jr^(\ + 2Vm;)(2p + ^fdW2(\ + pf){p + VM^(1 + p")) 

eQVo~2{-i + p?p 

Ml + 2^/M^)(2p + VW2{1 + P?){p + VW2{1 + p'))\T 



Oo^/eii-i + pfp 

The geometric scaling factor is given by iS. 8\) with 7 = {y/9i92)^^ ■ This design is 
also local E-optimal on the design space U = [s,t\ (0 < s < t), if the inequalities 
< s < p^^\/9i/62 and t > p^J 61/62 are satisfied. 

The local E-optimal design on the design space lA = [s,t\ for model h2. 2a\) is of 
the form ^3.9^) . if the inequalities s > p^^ \/ 61/62 and t > p\J 61/62 hold, of the 
form h3.1U\) . if the inequalities s < p^^ y/61/62 and t < p\J 61/62 are satisfied, 
and of the form Ii3.11\) . if the inequalities s > p^"^ \J 61/62 and t < p^/ 61/62 hold. 

5. Further discussion 

In this Section we discuss some practical aspects of the local optimal designs 
derived in the previous sections. In particular, we calculate the efficiency of a 
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design, which has recently been used in practice and investigate the efficiency 
of local optimal designs with respect to other optimality criteria. Throughout 
this paper the efficiency of a design ,^ is defined by eff$(^) = sup^ 
where $ denotes the particular optimality criterion under consideration and the 
optimal design maximizes 

Landete-Castillejos and Gallego (2000) used the inverse quadratic regression 
model to analyze data, which were obtained from lactating red deer hinds (Cervus 
elaphus). They concluded that inverse quadratic polynomials with parameteri- 
zation (|2.2ap can adequately describe the common lactation curves. The design 
space was given by the interval U = [1, 14], and the design used by these authors 
was a uniform design with support points (1, 2, 3, 4, 5, 6, 10, 14), which is denoted 
by throughout this section. The estimates for the parameters of model ()2.2ap 
are given by = 0.0002865, Oi = 0.0002117 and 62 = 0.0000301. Table 5.1 
shows the local optimal designs for the different optimality criteria considered 
in Section 3 and 4, where we used the point Xe = 21 for the calculation of the 
optimal extrapolation design. 



Criterion 


Optimal design 


D 


points : 


1 


3.4089 


14 


weights : 


1/3 


1/3 


1/3 


E 


points : 


1 


3.3561 


14 


weights : 


0.3972 


0.3914 


0.2114 


Di 


points : 


1 


3.3561 


14 


weights : 


0.1239 


0.2884 


0.5877 




points : 


1 


3.3561 


14 


Ce 


weights : 


0.0582 


0.1535 


0.7883 



Table 5.1. D-, E-, Di- and Ce-optimal designs for parametrization 112. 2a\) . 

The efficiencies of the different designs are shown in Table 5.2. We observe 
that the design of Landete-Castillejos and Gallego (2000) yields rather low effi- 
ciencies with respect to all optimality criteria, and the efficiency of the statistical 
analysis could have been improved by allocating observations according to local 
optimal design [see the first row in Table 5.2]. For example a confidence in- 
terval based on the local Di-optimal design would yield 66% shorter confidence 
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intervals for the parameter c as the design actually used by Landete-Castillejos 
and Gallego (2000). The advantages of the local optimal designs are also clearly 
visible for the other criteria. 





D 


E 


Di 


Ce 




69.92 


50.33 


45.85 


33.82 




100 


94.18 


75.28 


43.60 




93.96 


100 


51.89 


25.71 




74.63 


53.05 


100 


80.40 




51.23 


33.24 


85.73 


100 



Table 5.2. Efficiencies of local optimal designs and the uniform design for 
the inverse quadratic model (parameterization ^2.2a\) ) with respect to various 
alternative criteria (in percent). The design space is the intervallA = [1, 14], and 
the estimates of the parameters are given by 6q = 0.0002865, 9i = 0.0002117 and 
62 = 0.0000301. The local extrapolation optimal design is calculated for the point 
Xe = 21. 

Note that the data is usually used for several purposes, for example for 
discrimination between a linear and a quadratic inverse polynomial and for ex- 
trapolation using the identified model. Therefore it is important that an optimal 
design for a specific optimality criterion yields also reasonable efficiencies with 
respect to alternative criteria, which reflect other aspects of the statistical analy- 
sis. In Table 5.2 we also compare the efficiency of a given local optimal design 
with respect to the other optimality criteria. For example, the local D-optimal 
design has efficiencies 94.18%, 75.28% and 43.60% with respect to the E-, Di- 
and Ce-optimality criterion, respectively. Thus this design is rather efficient for 
the Di- and E'-optimality criterion, but less efficient for extrapolation. The situa- 
tion for the I?i -optimal design is similar, where the role of the Cg- and E'-criterion 
have to be interchanged. On the other hand the performance of the local E- and 
Ce-optimal design depends strongly on the underlying optimality criterion. The 
local ^'-optimal design yields only a satisfactory Z?-efficiency, but is less efficient 
with respect to the Cg- and Di-optimality criterion, while the local Cg-optimal 
design yields only a satisfactory Di-efficiency. 
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